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ABSTRACT 

We report on a 67 ks HETG observation of the optically brightest early 0-star, ^ Pup (04 f). Many 
resolved X-ray lines are seen in the spectra over a wavelength range of 5 to 25 A. Chandra has sufficient 
spectral resolution to study the velocity structure of isolated X-ray line profiles, and to distinguish the 
individual forbidden, intercombination, and resonance {fir) emission lines in several He-like ions even 
where the individual components are strongly Doppler broadened. In contrast with X-ray line profiles in 
other hot stars, C Pup shows blue-shifted and skewed line profiles, providing the clearest and most direct 
evidence that the X-ray sources are embedded in the stellar wind. The broader the line, the greater the 
blueward centroid shift tends to be. The N VII line at 24.78 A is a special case, showing a flat-topped 
profile. This indicates it is formed in regions beyond most of the wind attenuation. The sensitivity of 
the He-like ion fir lines to a strong UV radiation field is used to derive the radial distances at which 
lines of S XV, Si XIII, Mg XI, Ne IX, and O VII originate. The formation radii correspond well with 
continuum optical depth unity at the wavelength of each line complex, indicating that the X-ray line 
emission is distributed throughout the stellar wind. However, the S XV emission lines form deeper in 
the wind than expected from standard wind shock models. 

Subject headings: X-rays: stars — stars: individual (C Pup) 
outfiows — stars: mass-loss — line: profiles 



stars: early-type — stars: winds. 



other OB stars. Wind shock models developed by Lucy 
& White (1980), Feldmeier et al. (1997) and others, con- 
sistently failed to predict the high levels of X-ray emis- 
sion observed in the brightest O stars like ^ Pup, lead- 
ing to the suggestion that perturbations somehow form 
and propagate up from the photosphere into the wind and 
drive stronger shocks (Feldmeier 1995; Feldmeier, Puis, 
& Pauldrach 1997). Broad-band X-ray observations of 
C Pup (Corcoran et al. 1993; Hillier et al. 1993) indi- 
cate that some wind attenuation is affecting the soft X-ray 
flux. However with the advent of Chandra and XMM we 
can apply, for the first time, diagnostic emission line ratios 
and measure line profiles to determine the locations and 
Doppler velocities of X-ray sources in the stellar wind of 
CPup. 

2. OBSERVATIONS AND CONSTRAINTS FROM 
HELIUM-LIKE fir LINE RATIOS 

We obtained a 67 ks Chandra High Energy Transmis- 
sion Grating Spectrometer (HETGS) observation of C Pup 
from 2000 Mar 28, 13'' 31™ UT to Mar 29, 09'' 12". 
The standard pipeline tools were used to reprocess the 
data with the most recent calibration files available. Line 
emission is clearly evident in both the high-energy grating 
(HEG) and the medium-energy grating (MEG). The com- 
bined ibl'** order, background-subtracted, HETGS spectra 
are shown in Figure 1. Triads of He-like ions known as 
the fir (forbidden, intercombination, and resonance) lines 
(S XV, Si XIII, Mg XI, Ne IX, and O VII) are seen, in ad- 
dition to isolated Lya emission lines and numerous L-shell 
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1. INTRODUCTION 

The 04f star C, Puppis has for decades been at the cut- 
ting edge of research regarding early type stars because 
of its optical and UV brightness. In addition, because of 
the low interstellar attenuation, C, Pup has been the prime 
X-ray target to study soft X-ray emission from O stars. 
Researchers have looked to this source to settle outstand- 
ing controversies about the physical location, quantity, and 
nature of the hot X-ray emitting plasma on OB stars. In 
this paper we present the highest resolution X-ray spec- 
trum of C Pup ever measured and explore the physical 
mechanism of hot star X-ray production. 

Prior to the discovery of X-ray emission from O stars 
(Harnden et al. 1979; Seward et al. 1979), Cassinelli 
& Olson (1979) postulated thin base coronal zones as a 
source for X-rays to explain the observed UV superion- 
ization (Lamers & Morton 1976) by way of the Auger ef- 
fect. The coronal models (Cassinelli & Olson 1979; Wal- 
dron 1984) predicted an X-ray absorption at the oxygen 
K-shell edge larger than that observed with Einstein SSS 
(Cassinehi & Swank 1983) and BBXRT (Corcoran et al. 
1993). These observations, as well as detailed modeling of 
the superionization profiles (MacFarlane et al. 1993) and 
lack of detection of the iron "green line" (Baade & Lucy 
1987) in C, Pup caused the coronal model to fall out of 
favor. Increasingly, a consensus formed around a wind- 
shock picture in which a series of shocks, perhaps related 
to the line-force instability (Lucy 1982; Owocki, Castor, 
& Rybicki 1988), causes hot. X-ray emitting gas to be dis- 
tributed throughout the dense stellar wind of C Pup and 
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lines of iron, especially Fe XVII. 
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Fig. 1. — Co-added 1"* order, background subtracted, Chandra 
HETGS HEG (top panel) and MEG (bottom panels) count spectra 
of C Pup. The ions responsible for the strongest line emissions are 
identified. The bin size is 0.01 A. 

Traditionally, He- like ion f /i line ratios have been used 
to derive electron densities of X-ray line emitting regions 
since the populations of the 2^P levels are controlled by 
collisional excitations from the 2''S levels (Gabriel & Jor- 
dan 1969). However, when a strong external UV source 
is present, the excitation 2'^S — > 2'^P is predominantly ra- 
diative and this means that the f/i ratio is no longer a 
density diagnostic, but rather it can be used to determine 
the strength of the UV radiation field (Blumenthal et al. 
1972). The UV wavelengths associated with the radia- 
tive excitation are in the range from ~ 650 to 1650 A, 
where O star UV fluxes are very strong. For ( Pup's pho- 
tosphcric flux wo use a Kurucz (1993) model, assmning 
reff= 42000K. Since the 2^8 wavelengths are not 

on the Wien side of the spectrum, the model atmosphere 
fluxes arc relatively well determined. 

Waldron & Cassinelli (2001) demonstrated how the ob- 
served f/i line ratios can bo used to derive the line for- 
mation radii from the geometric dilution of the photo- 
spheric radiation field in ( Ori (09.5Ia). Kahn et al. (2001) 
used a similar approach to derive line formation radii 
from their XMM RGS spectrum of C, Pup. Our Chandra 
HETG spectrum has better resolution at essentially all He- 
like ion wavelengths and greater sensitivity to the short- 
wavelength lines (e.g., Si and S) than the XMM spectrum. 
This high resolution is very important for O stars since 
their X-ray lines show strong broadening. The observed 
He- like f/i line ratios are listed in Table 1. 



Figure 2a shows the radii derived from our measure- 
ments of the He-like f/i ratios. The radii of line formation 
are seen to vary from ion to ion, with lines from high nu- 
clear charge ions (Si, S) forming at small radii, and lower 
nuclear charge ion lines (O, Ne, Mg) forming at larger 
radii. One might expect the fir line complexes for each 
element to be formed over a range of radii and not just in a 
thin shell. However, since the line emission is proportional 
to n^, the observed emission in each line complex should 
be (loniinatcd by the densest regions (i.e. the smallest ra- 
dius) from which the line radiation can escape. For Q Ori 
the fir formation locations correlates with the associated 
radial optical depth unity radius, Ri. The optical depth is 
measured through the wind, using the continuum opacity 
at the wavelength of each triad of fir lines. Any line ra- 
diation that originates much below R\ could not be seen 
by an external observer. The wavelength dependence of 
Ri for C Pup is shown in Figure 2b using a /3-law with 
(3 = 0.75 as the velocity law, the wind absorption cross 
sections in Waldron et al. (1998), and the stellar param- 
eters of Lamers & Leithcrer (1993) (i?* = 16 Rq, M= 
2.4 X 10"*^ Meyr-i, and v-y,= 2200 kms^^). 
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Fig. 2. — The top panel shows the dependence of O VII, Ne IX, 
Mg XI, Si XIII and S XV f/i ratios on radius due to the geometric 
dilution of the UV radiation field. The length of the darkened line 
for each ion extends over the range of the uncertainty in the f/i line 
ratio. The bottom panel gives the radius of wind continuum opti- 
cal depth unity as a function of wavelength. The dark horizontal 
lines correspond to the range in radii indicated in the upper panel. 
For comparison, we also show the optical depth unity radius assum- 
ing neutral interstellar medium opacity of Morrison & McCammon 
(1983). 

An inspection of Figure 2 leads to the following conclu- 
sions: First, the He-like line formation radii are all consis- 
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Table 1: fir ratio measurements 



SXV 


Si XIII 


Mg XI 


Ne IX 


VII 


f/i 0.61 ± 0.48 


1.03 ±0.14 


0.29 ±0.18 


< 0.1 


< 0.1 


r/i 0.78 ± 0.54 


2.00 ±0.25 


1.53 ±0.13 


1.44 ±0.18 


1.30± 0.6 


J^tot" 0.81 


3.06 

in— 13 


3.26 

-2 „-l 


5.77 


3.56 



"Total flux of the three lines in 10 ^'^ erg cm 



tent with their associated Ri values, as expected from the 
above argument. Second, the f/i ratio of S XV indicates 
that hot plasma exists below 1.2 R* where the nominal 
wind velocity is < 600 kms~^. Most radiative instability 
simulations do not show strong shocks this close to the star 
(Feldmeier 1995, Owocki, Castor, & Rybicki 1988). How- 
ever, using a velocity jump even half as large as the local 
wind velocity (Av = 300km s"-'^) to calculate the shock 
temperature from T = 1.4 x lO^if (Ai;/100kms~^)^, we 
obtain a post-shock temperature of 1 MK; too low to pro- 
duce S XV. Thus, anomalously strong shock jumps are re- 
quired deep in the wind. A similar conclusion was reached 
by Waldron & CassincUi (2001), but for ( Ori the highest 
ion stage observed was Si XIII instead of the S XV seen 
here. 

3. EMISSION LINE PROFILES 

The most notable feature observed in our HETG spec- 
trum of ( Pup is the clear presence of blue-shifted X-ray 
line centroids in all strong lines. This is in contrast with 
previous Chandra O star observations (Schulz et al. 2000; 
Waldron & Cassinelli 2001) where the lines were broad, 
symmetric and un-shifted. Although our line widths are 
comparable to other 0-star observations, with HWHM of 
~ 1000 kms^^, our observation is the first Chandra detec- 
tion of blue-shifted and blueward skewed X-ray line profiles 
in an O star. Figure 3 displays six ( Pup X-ray lines ar- 
ranged in order of wavelength. The results for the line fits 
are given in Table 2. There is an interesting progression 
that holds for these lines: the centroid shifts are generally 
larger for the broader lines, indicating a connection be- 
tween radial location, wind absorption, and the Doppler 
broadening of the line emission. In Table 4 of their XMM 
study of C Pup, Kahn et al. (2001) list blueshifts of the 
Lyman-a lines of Ne X and O VIII. The blueshifts agree 
with ours to within 2 sigma, and some of the discrep- 
ancy can be attributed to the absohitc wavelength cali- 
bration uncertainties of the instruments 100 km 
for the HETG (Chandra X-ray Center (CXC) 2001), see 
den Herder et al. 2001 for the RGS). The XMM results 
show the same trend between shift and broadening. The 
N VII line is an exception to the trend in Table 2 because 
of its small centroid shift. It is different morphologically 
from the other lines (also noted by Kahn et al. (2001)). It 
has a roughly flat-topped profile so we use a box fit instead 
of a Gaussian to estimate the values for Table 2. 

There is a blueward skewness observed in many of the 
( Pup lines which is not attributable to the HETGS line 
response function (CXC 2001). The blue sides of these 
lines tend to be steeper, while the red sides have a shal- 
lower slope. In general, the shape of these lines indicates 
wind broadening combined with attenuation as would arise 
from line formation in the marginally thick region near 
Ri. MacFarlane et al. (1991), Ignace (2001), and Owocki 



& Cohen (2001) demonstrate that emission line profiles 
can provide important information about the spatial dis- 
tribution of hot plasma within a stellar wind based both 
on the velocity-dependence of the intrinsic emission and 
also on the continuum attenuation across the line. As- 
suming isotropic emission, the red-shifted emission from 
the back side of the wind is suppressed by the continuum 
opacity along the line of sight, while the blue-shifted emis- 
sion from the near side is less so. This leads to skewed, 
almost triangular-shaped lines in an optically thick wind. 
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Fig. 3. — Six X-ray line profiles from the Pup spectrum (solid) 
with fits (dashed). The vertical dotted lines indicate the rest wave- 
lengths for these transitions. The first 5 lines were fit with Gaussians 
while the N VII line was fit with a box function because of its fiat- 
topped shape. All functions shown have been convolved with the 
instrumental response. These functions are simultaneous fits to the 
MEG+1 and MEG-1 order spectra but are shown here superimposed 
on the co-added MEG first order spectrum for the sake of brevity. 
The data have been binned to 0.01 A. The fit parameters are given 
in Table 2. 

In the C Ori HETG spectrum, Waldron & Cassinelli 
(2001) found that an optically thin model was needed to 
provide a reasonable fit because the lines were symmetric 
and un-shifted. This fit required that the assumed mass 
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Table 2: Line Profile Widths and Centroid Shifts 



Ion 


Arest 


•^line 


Source 


Centroid 










HWHM 


Shift 






/ X \ 

(A) 




(km s ) 


(km/sj 


[R*) 


Mg XII 


8.42 


0.30 


610 ±220 


-440 ± 130 


2.8 


NeX 


12.13 


3.33 


780 ± 60 


-570 ± 50 


4.5 


Fe XVII 


15.01 


4.74 


800 ± 70 


-500 ±60 


6.0 


Fe XVII 


16.78 


2.24 


850 ± 130 


-730 ± 100 


7.0 


VIII 


18.97 


3.05 


990 ± 201 


-670 ± 140 


4.0 


N VII 


24.78 


5.02 


1570 ± 150 


-120 ±150 


4.0 



"Line flux in 10 ergs cm ^ s 



loss rate for ^ Ori be greatly reduced. For C. Pup we have 
the odd situation that the blueward shift and asymmet- 
ric line profiles agree with what has long been expected 
for hot stars, but this is the first time that it has been 
so clearly observed in a Chandra spectrum for any star. 
We find that a fit can be made to the line profiles without 
making a significant reduction in the assumed mass loss 
rate in contrast with the ( Ori results. 

The N VII lino profile is different from those of the other 
ions as can be seen in Figure 3 and Table 2. The Doppler 
half- width of this line is much greater than any of the other 
X-ray lines. This velocity of 1570 kms^^ can be compared 
with the terminal speed of 2200 kms^^ determined from 
UV wind lines. A flat-topped line profile such as this is 
expected from a fast moving shell source suffering little or 
no wind attenuation (MacFarlane et al. 1991). Therefore, 
in contrast with the other observed lines, the NVII profile 
is probably formed well above its continuum optical depth 
unity radius. 

4. DISCUSSION AND SUMMARY 

The Chandra HETG emission line spectrum of ( Pup, 
with its unprecedented resolution (in excess of ^ = 1000 
for some lines) and sensitivity over a wide range of wave- 
lengths, allows us to draw quantitative conclusions about 
the nature of the X-ray source on this prototypical hot 



star. Most importantly, we find that the wind plays an 
important, observable role in determining the X-ray line 
profiles of ( Pup, in contrast to results for other stars. 

We can begin to constrain the physical processes which 
lead to the production of the hot, X-ray emitting plasma 
on this star. Taking all the He-like ions from oxygen 
through sulfur into account, the correlation between the 
formation radius of the fir lines and the radius at opti- 
cal depth unity, -Ri, we find evidence for a spatially dis- 
tributed source of X-rays throughout the expanding wind. 

There are two lines of particular interest in the develop- 
ment of future shock models: S XV, and N VII. The shock 
jump required to produce S XV appears inconsistent with 
the expected local wind conditions, requiring a post-shock 
flow velocity of approximately zero. Although this is not 
impossible for a shock model to reproduce, it places strong 
constraints on future modeling efforts. The N VII Lya line 
line may have a shape distinct from the other lines in our 
spectrum because it is formed in the outer regions of the 
wind where there is no further acceleration of the X-ray 
source region and little overlying wind attenuation. This 
Chandra HETG data set provides the most detailed and 
complete picture to date of an extended distribution of hot 
plasma embedded within a strong, optically thick stellar 
wind. 
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